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Abstract 

The quantum Zakharov system in three-spatial dimensions and an associated Lagrangian descrip- 
tion, as well as its basic conservation laws are derived. In the adiabatic and semiclassical case, 
the quantum Zakharov system reduces to a quantum modified vector nonlinear Schrodinger (NLS) 
equation for the envelope electric field. The Lagrangian structure for the resulting vector NLS 
equation is used to investigate the time-dependence of the Gaussian shaped localized solutions, 
via the Rayleigh-Ritz variational method. The formal classical limit is considered in detail. The 
quantum corrections are shown to prevent the collapse of localized Langmuir envelope fields, in 
both two and three-spatial dimensions. Moreover, the quantum terms can produce an oscillatory 
behavior of the width of the approximate Gaussian solutions. The variational method is shown to 
preserve the essential conservation laws of the quantum modified vector NLS equation. 
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I. INTRODUCTION 



The Zakharov system [jj, describing the couphng between Langmuir and ion-acoustic 
waves, is one of the basic plasma models, see Ref. [2I, S] for reviews. Recently 4], a quantum 
modified Zakharov system was derived, by means of the quantum plasma hydrodynamic 
model In this context, enhancement of the quantum effects was then shown e. g. 

to suppress the four- wave decay instability. Subsequently j^, a kinetic treatment of the 
quantum Zakharov system has shown that the modulational instability growth rate can be 
increased in comparison to the classical case, for partially coherent Langmuir wave electric 
fields. Also [9], a variational formalism was obtained and used to study the radiation of 
localized structures described by the quantum Zakharov system. Bell shaped electric field 
envelopes of electron plasma oscillations in dense quantum plasmas obeying Fermi statistics 
were analyzed in Ref. 1^. More mathematically-oriented works on the quantum Zakharov 
equations concern its Lie symmetry group and the derivation of exact solutions 
Finally, there is evidence of hyper chaos in the reduced temporal dynamics arising from the 



quantum Zakharov equations 



15|- 



All these paper refer to quantum Zakharov equations in one-spatial-dimension only. In 
the present work, we extend the quantum Zakharov system to fully three-dimensional space, 
allowing also for the magnetic field perturbation. In the classical case, both heuristic argu- 
ments and numerical simulations indicate that the ponderomotive force can produce finite- 
time collapse of Langmuir wave packets in two- or three-dimensions j^, [16, 17 1. This is in 



contrast to the one-dimensional case, whose solutions are smooth for all time. A dynamic 
rescaling method was used for the time-evolution of electrostatic self-similar and asymp- 
totically self-similar solutions in two- and three-dimensions, respectively 18|. Allowing for 



transverse fields shows that singular solutions of the resulting vector Zakharov equations 
are weakly anisotropic, for a large class of initial conditions [l9|. The electrostatic nonlin- 
ear collapse of Langmuir wave packets in the ionospheric and laboratory plasmas has been 



observed 
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2l| . Also, the collapse of Langmuir wave packets in beam plasma experiments 



verifies the basic concepts of strong Langmuir turbulence, as introduced by Zakharov |22| . 
The analysis of the coupled longitudinal and transverse modes in the classical strong Lang- 



muir turbulence has been less studied 



23j-25|], as well as the intrinsically magnetized case 



261], which can lead to upper- hybrid wave collapse [27|]. Finally, Zakharov-like equations 
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have been proposed for the electromagnetic wave collapse in a radiation background 28|. 

It is expected that the ponderomotive force causing the collapse of localized solutions 
in two- or three-space dimensions could be weakened by the inclusion of quantum effects, 
making the dynamics less violent. This conjecture is checked after establishing the quan- 
tum Zakharov system in higher-dimensional space and using its variational structure in 
association with a (Rayleigh-Ritz) trial function method. 

The manuscript is organized in the following fashion. In Section 2, the quantum Zakharov 
system in three-spatial-dimensions is derived by means of the usual two-time scale method 
applied to the fully 3D quantum hydrodynamic model. In Section 3, the 3D quantum Za- 
kharov system is shown to be described by a Lagrangian formalism. The basic conservation 
laws are then also derived. When the density fluctuations are so slow in time so that an 
adiabatic approximation is possible, and treating the quantum term of the low-frequency 
equation as a perturbation, a quantum modified vector nonlinear Schrodinger equation for 
the envelope electric field is obtained. In Section 4, the variational structure is used to ana- 
lyze the temporal dynamics of localized (Gaussian) solutions of this quantum NLS equation, 
through the Rayleigh-Ritz method, in two-spatial-dimensions. Section 5 follows the same 
strategy, extended to fully 3D space. Special attention is paid to the comparison between 
the classical and quantum cases, with considerable qualitative and quantitative differences. 
Section 6 contains conclusions. 



II. QUANTUM ZAKHAROV EQUATIONS IN 3 + 1 DIMENSIONS 

The starting point for the derivation of the electromagnetic quantum Zakharov equations 
is the quantum hydrodynamic model for an electron- ion plasma. Equations (20)- (28) of Ref. 
[3]. For the electron fluid pressure pe, consider the equation of state for spin 1/2 particles 
at zero temperature, 

o 2 5/3 

= (1) 

where me is the electron mass, vp^ is the Fermi electron thermal speed, ng is the electron 
number density and is the equilibrium particle number density both for electron and 
ions. The pressure and quantum effects (due to their larger mass) are neglected for the ions. 
Also due to the larger ion mass, it is possible to introduce a two-time scale decomposition, 
= '^o + <^^s + <^%) ''^i = "".o + ^ris-, Ue = (Ju^ + (5uj, Uj = (Ju^, E = (JE^ + (5Ej, B = 5Bj, 



where the subscripts s and / refer to slowly and rapidly changing quantities, respectively. 
Also, Ug is the electron fluid velocity, rii the ion number density, Uj the ion fluid velocity, E 
the electric fleld, and B the magnetic fleld. Notice that it is assumed that there is no slow 
contribution to the magnetic fleld, a restriction which allows to get B = (mg/e) V x 6uf (see 
Equation (2.21) of Ref. ^]), where — e is the electron charge. Including a slow contribution 
to the magnetic fleld could be an important improvement, but this is outside the scope of 
the present work. 

Following the usual approximations [3|, the quantum corrected 3D Zakharov equations 
read 

2iu;pe— - V X (V X E) + vl^V{V ■ E) = 



2 



— tij^ EH V 

Uq Aral 



V2(V-E) , (2) 



^ - WHus - -^V^dEp) + W'5ns = . (3) 

Here E is the slowly varying envelope electric fleld deflned via 

E/ = i (E e"'"^^^* + E* e'''^^*) , (4) 

where uope is the electron plasma frequency. Also, in Eqs. (I2H3]) c is the speed of light in 
vacuum, h the scaled Planck constant, Eq the vacuum permittivity and mi the ion mass. 
In addition, = KBTpf./mi , where KbTfc = ^ev'p^- Therefore, Cg is a Fermi ion-acoustic 
speed, with the Fermi temperature replacing the thermal temperature for the electrons. 

In comparison to the classical Zakharov system (see Eqs. (2.48a)-(2.48b) of Ref. j^), 
there is the inclusion of the extra dispersive terms proportional to It? in Eqs. (l2])-(l3l). 
Other quantum difference is the presence of the Fermi speed instead of the thermal speed 
in the last term at the left hand side of Eq. ([2]). From the qualitative point of view, the 
terms proportional to are responsible for extra dispersion which can avoid collapsing of 
Langmuir envelopes, at least in principle. This possibility is investigated in Sections 4 and 
5. Finally, notice the non trivial form of the fourth order derivative term in Eq. ([2]). It is 
not simply proportional to V^E as could be wrongly guessed from the quantum Zakharov 
equations in 1 + 1 dimensions, where there is a ~ d'^'Ei/dx'^ contribution ^. 
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It is useful to consider the rescaling 



t = 2fl Upet , (5) 



e 



Sris ^ eE 
n = , t 



where = rrie/mi. Then, dropping the bars in r, t, we obtain 

.d£ c2 



- X (V X ^) + V(V ■ 8) 

ot vi 



= nS + TV [V^(V-^)] , (6) 



where 



V'n-V\\S\')+TV^n = 0, (7) 



hujpf, ^ 



TTLi \KBTFeJ 

is a non-dimensional parameter associated with the quantum effects. Usually, it is an ex- 
tremely small quantity, but it is nevertheless interesting to retain the ~ F terms, specially 
for the collapse scenarios. The reason is not only due to a general theoretical motivation, 
but also because from some simple estimates one concludes that these terms become of the 
same order as some of other terms in Eqs. ([S])^© provided that the characteristic length / 

— 1/3 

for the spatial derivatives becomes as small as the mean inter-particle distance, / ~ tzq . 
Of course, the Zakharov equations are not able to describe the late stages of the collapse, 
since they do not include dissipation, which is unavoidable for short scales. But even Landau 
damping would be irrelevant for a zero-temperature Fermi plasma, where the main influence 
comes from the Pauli pressure. In the left-hand side of Eq. ([6]), the V(V-£^) term is retained 
because the ~ (? I'^\e transverse term disappears in the electrostatic approximation. 

In the adiabatic limit, neglecting d'^n/dt^ in Eq. ([7]) and under appropriated boundary 
conditions, it follows that 

n = -l^p + rV^n, (9) 



When F 7^ 0, it is not easy to directly express n as a function of \£\ as in the classical 
case. Therefore, the adiabatic limit is not enough to derive a vector nonlinear Schrodinger 
equation, due to the coupling in Eq. ([9]). 
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III. LAGRANGIAN STRUCTURE AND CONSERVATION LAWS 



The quantum Zakharov equations ([6])-([7]) can be described by the Lagrangian density 



2\ dt dt 



|v X |v-fp-r|v(v-f)p 



rlVnl 



iVal 



(10) 



where n, the auxihary function a and the components of S,S* are regarded as independent 
fields. Remark: for the particular form (fTOj) and for a generic field ip, one computes the 
functional derivative as 



6C dC d dC 



d dC 

+ 



dC 



5ip dip dri dip/ Or i dt dip/dt dridrj d'^ip/dridrj 

using the summation convention and where are cartesian components. 
Taking the functional derivatives with respect to n and a, we have 

— = n + \c\ — rV n, 
dt 



(11) 



and 



dn 



(12) 



(13) 



respectively. Eliminating a from Eqs. (1121) and fll3p we obtain the low frequency equation. 
In addition, the functional derivatives with respect to S* and S produce the high-frequency 



equation and its complex conjugate. The present 



formulation of the classical Zakharov equations 29[ |. 



brmalism is inspired by the Lagrangian 



The quantum Zakharov equations admit as exact conserved quantities the "number of 
plasmons" of the Langmuir field, 

N= f \£\^dr, (14) 



the linear momentum (with components Pi, i = x,y, z) 



l2 



d_si 

dvi 



d£^ 

dvi 



— n 



da 
dvi 



dr 



(15) 



and the Hamiltonian, 

n -- 



+ — |V X + |v-^p + r|v(v-^)p 



+ 



•'Fe 



n2 + r|Vnp+ iVaP 



dr , 



(16) 



Furthermore, there is also a preserved angular momenta functional, but it is not relevant 
in the present work. These four conserved quantities can be associated, through Noether's 
theorem, to the invariance of the action under gauge transformation, time translation, space 
translation and rotations, respectively. The conservation laws can be used e. g. to test the 
accuracy of numerical procedures. Also, observe that equations ([7]) and for the adiabatic 
limit are described by the same Lagrangian density (fTOj) . In this approximation, it suffices 
to set a = 0. 

In addition to the adiabatic limit, Eq. (Q can be further approximated to 



n 



-\S\ 



(17) 



assuming that the quantum term is a perturbation. In this way and using Eq. 
quantum modified vector nonlinear Schrodinger equation is derived 



©, a 



tj- + V(V-^)-— Vx (Vx^) 

= rv [v'(v-f)] -r^v'd^i^ 



(18) 



The appropriate Lagrangian density Cad,sc for the semiclassical equation f[T5]) is given by 



ad,sc 



-(£*■ — 

2 y"- dt 



r|v(v-^) 



dS* 
'dt 



^|V X S\ 



|V-£| 



r 

2 



V[|^| 



(19) 



where the independent fields are taken as £ and £* components. 

The expression for the number of plasmons in Eq. f fT4l) remains valid as a constant 
of motion in the joint adiabatic and semiclassical limit, as well as the momentum P in Eq. 
( TTSl) with a = 0. Finally, the Hamiltonian 

^2 



ad.sc 



^|V X S\' + \V -81' + T\V{V -8) 



2 ' ' 2 



Vil^l 



dr 



(20) 



is also a conserved quantity. 

In the following, the infiuence of the quantum terms in the right-hand side of Eq. flTS]) 
are investigated, assuming adiabatic conditions for collapsing quantum Langmuir envelopes. 
Other scenarios for collapse, like the supersonic one 
be investigated in the future. 
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19| , could also be relevant and shall 



IV. VARIATIONAL SOLUTION IN TWO DIMENSIONS 



Consider the adiabatic semiclassical system defined by Eq. (fT8l) . We refer to locahzed 
solution for this vector NLS equation as (quantum) "Langmuir wave packets" , or envelopes. 



As discussed in detail in 



29| in the purely classical case, Langmuir wave packets will become 



singular in a finite time, provided the energy is not bounded from below. Of course, explicit 
analytic Langmuir envelopes are difficult to derive. A fruitful approach is to make use of the 
Lagrangian structure for deriving approximate solutions. This approach has been pursued 



m 



30l | for the classical and in for the quantum Zakharov system. Both studies considered 



the internal vibrations of Langmuir envelopes in one-spatial-dimension. Presently, we shall 
apply the time-dependent Rayleigh-Ritz method for the higher-dimensional cases. A priori, 
it is expected that the quantum corrections would inhibit the collapse of localized solutions, 
in view of wave-packet spreading. To check this conjecture, and to have more definite 
information on the infiuence of the quantum terms, first we consider the following Ansatz, 

— ] - exp f-^j exp(z(e + V)) (cos0,sin0,O), (21) 

which is appropriate for two-spatial-dimensions. Here cr, k, G and are real functions of 
time, and p = x'^ + y"^. The normalization condition f|T^ is automatically satisfied (in 
2D the spatial integrations reduce to integrations on the plane). Other localized forms, 
involving a seek type dependence, conld have been also proposed. Here a Gan=.ian 

form was suggested mainly for the sake of simplicity |31|. Notice that the envelope electric 
field ( 12T|) is not necessarily electrostatic: it can carry a transverse (V x £^ 7^ 0) component. 

The free functions in Eq. f l2T|) should be determined by extremization of the action 
functional associated with the Lagrangian density (fT9|) . A straightforward calculation gives 



^ad,scdxdy = -N 



■ 2; 1 / c2 N' 



+ SrP + lOTkV + (22) 

\ Zn J J 

where only the main quantum contributions are retained. Now L2 is the Lagrangian for a 
mechanical system, after the spatial form of the envelope electric field was defined in advance 
via Eq. fl2T|) . Of special interest is the behavior of the dispersion a. For a collapsing solution 
one could expect that a goes to zero in a finite time. The phase G and the chirp function 
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k should be regarded as auxiliary fields. Notice that L2 is not dependent on the angle </>, 
which remains arbitrary as far as the variational method is concerned. 
Applying the functional derivative of L2 with respect to 6, we obtain 

^ = - iV = 0, (23) 

so that the variational solution preserves the number of plasmons, as expected. The remain- 
ing Euler-Lagrange equations are 



ok V 



^ aa = — a^A; + SrA; + 32raV , (24) 

"'Pe 



SL2 : 2c2 „ 1 / c2 \ 1 . 



5a " 2 V^Fe 27ry a^ 

The exact solution of the nonlinear system (!24] - [25|) is difficult to obtain, but at least the 
dynamics was reduced to ordinary differential equations. 

It is instructive to analyze the purely classical (F = 0) case first. This is specially true, 
since to our knowledge the Rayleigh-Ritz method was not applied to the vector NLS equation 
(lisp , even for classical systems. The reason can be due to the calculational complexity 
induced by the transverse term. When F = 0, Eq. (p^ gives k = Vpg&/2c^a. Inserting this 
in Eq. fl25l) we have 

where the pseudo-potential V2C is 

2"^ \v^,~2^)V^ 

From Eq. (1271) it is evident that the repulsive character of the pseudo-potential will be 
converted into an attractive one, whenever the number of plasmons exceeds a threshold, 

Ar>^, (28) 

a condition for Langmuir wave packet collapse in the classical two-dimensional case. The 
interpretation of the result is as follows. When the number of plasmons satisfy Eq. (1251) . 
the refractive ~ \£\^ term dominates over the dispersive terms in the Lagrangian density 
( IT9|) . producing a singularity in a finite time. Finally, notice the ballistic motion when 
N = 27ic^/vp^, which can also lead to singularity. 



Further insight follows after evaluating the energy integral (120|) with the Ansatz (12T!) . 
which gives, after eliminating k, 

,2 



K 



ad,sc,2c 



(r^o). 



(29) 



Of course, this energy first integral could be obtained directly from Eq. fl26l) . However, the 
plausibility of the variational solution is reinforced, since Eq. ( l29l) shows that it preserves 
the exact constant of motion Tiad^sc- In addition, in the attractive (collapsing) case the 
energy (l29l) is not bounded from bellow. 

In the quantum (F 7^ 0) case, Eq. (p^ becomes a cubic equation in fc, whose exact 
solution is too cumbersome to be of practical use. It is better to proceed by successive 
approximations, taking into account that the quantum and electromagnetic terms are small. 
In this way, one arrives at 



a 



' da ' 



where the pseudo-potential V2 is 



2?;2 



N 
2^ 



Fc2 



-'Fe 



1 + 



N\ 1 



- 



27ry 



(30) 



(31) 



Now, even if the threshold ( |28l) is exceeded, the repulsive ~ cr~^ quantum term in V2 will 
prevent singularities. This adds quantum diffraction as another physical mechanism, besides 
dissipation and Landau damping, so that collapsing Langmuir wave packets are avoided in 
vector NLS equation. Also, similar to Eq. (!29|) . it can be shown that the approximate 
dynamics preserves the energy integral, even in the quantum case. Indeed, calculating from 
Eq. fl2UI) and the variational solution gives "hiad sc 

^2 



n 



ad,sc,2 



a 



+ V2 



:f>o). 



(32) 



From Eq. (1501) . obviously 'Had,sc,2 = 0. 

It should be noticed that oscillations of purely quantum nature are obtained when the 
number of plasmons exceeds the threshold fl28|) . Indeed, in this case the pseudo-potential 
V2 in Eq. (I3T]) assumes a potential well form as shown in Figure 1, which clearly admits 
oscillations around a minimum a = am- Here, 

1/2 



T{1 + N/2tt) 
[N/271 -cyvl^ 



(33) 
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Also, the minimum value of V2 is 



c2 (iV/27r - cV4e)' ^ (N c2^' 



V2{(ym) = - .^p 2 1 , Ar/o ~ > ~THF 2~ ' (^'^) 

loi l+AZ/ZTT loi yZTT "W^e/ 

the last inequality follows since Eq. ( l28l) is assumed. Therefore, a deepest potential well 
is obtained when is increasing. Also, for too large quantum effects the trapping of the 
localized electric field in this potential well would be difficult, since V2{<7m) 0_ as F 
increases. This is due to the dispersive nature of the quantum corrections. 

The frequency u of the small amplitude oscillations is derived linearizing Eq. ( !30|) around 
the equilibrium point fl5^ . Restoring physical coordinates via Eq. this frequency is 
calculated as 

c f KBTpe V (Ar/27r - d'/vhf/^ ^ 

<^f^Vf^-:^V^.. (35) 



^/2c V ^^pe J V27r 



V 



Fe 



To conclude, the variational solution suggests that the extra dispersion arising from 
the quantum terms would inhibit the collapse of Langmuir wave packets in two-spatial- 
dimensions. Moreover, for sufficient electric field energy (which is proportional to A^), in- 
stead of collapse there will be oscillations of the width of the localized solution, due to the 
competition between the classical refraction and the quantum diffraction. The frequency of 
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linear oscillations is then given by Eq. (l35l) . The emergence of a pulsating Langmuir envelope 
is a qualitatively new phenomena, which could be tested quantitatively in experiments. 



V. VARIATIONAL SOLUTION IN THREE-DIMENSIONS 



It is worth to study the dynamics of localized solutions for the vector NLS equation f|T8|) 
in fully three-dimensional space. For this purpose, we consider the Gaussian form 



N 



1/2 

>3 / 



2^2 



fcos d) sin 9, sin 6 sin 9, cos ( 



(36) 



where cr, k, 6, 9 and are real functions of time and r = + y"^ + z^, applying the 
Rayleigh-Ritz method just like in the last Section. The normalization condition (IT^ is 
automatically satisfied with Eq. fl36l) . which, occasionally, can also support a transverse 
(V X ^ ^ 0) part. 

Proceeding as before, the Lagrangian 

N 



+ 10 r/t^ + 20 TA; V 



sr 



3TN 



(37) 



4V2 7r3/V5 

is derived. In comparison to the reduced 2D-Lagrangian in Eq. ( 122|) . there are different 
numerical factors as well as qualitative changes due to higher-order nonlinearities. Also, the 
angular variables 9 and don't appear in L3. 

The main remaining task is to analyze the dynamics of the width a as a function of time. 
This is achieved from the Euler-Lagrange equations for the action functional associated to 
L3. As before, 6L3/6Q = gives = 0, a consistency test satisfied by the variational 
solution. The other functional derivatives yield 



6L3 
5k 

6a 








ak 



Ak 

y 

1 

3 



2c2 



2_4\ 



a' + sr (1 + Ak'a 
2 o ^2 



(38) 



i , o 2 c 



3A^ 



- sorra 



-"Fe 



4„3 



4y2 7r3/2 



sr 



IhVN 



4y2 7r3/2a6J 



In the formal classical limit (r = 0), and using Eq. (1381) to eliminate k, we obtain 



a 



da 



(39) 



(40) 
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where now the pseudo-potential V^c is 

^ r_8^ 

4e U^^Fe^' 9V2 7r3/2a3 



3c 



2iV 



The form (l4Tll shows a generic singular behavior, since the attractive 



~ 0" 



(41) 



term will 



dominate for sufficiently small a, irrespective of the value of A^. Hence, in fully three- 
dimensional space there is more "room" for a collapsing dynamics. Figure 2 shows the 
qualitative form of V^c, attaining a maximum at o" = ctm, where 



3vlN 



(42) 



8v^7r3/2c2 ■ 

By Eq. f l39|) and using successive approximations in the parameter F to eliminate k via 
Eq. (l38!l . we obtain 

(43) 



da 



where 



8c^ 
3^ 



5F 

+ + 



FA^ 



f44) 

.34,cr2 12v^7r3/2a3 Ua^ 4^2713/2^5] ' ^ ^ 

The quantum terms are repulsive and prevent collapse, since they dominate for sufficiently 
small a. Moreover, when F 7^ an oscillatory behavior is possible, provided a certain 
condition, to be explained in the following, is meet. 

To examine the possibility of oscillations, consider V^{a) = 0, the equation for the critical 
points of V3. Under the rescaling s = cr/aM, where (Tm (defined in Eq. (H2l) ) is the maximum 
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cr 



cr 



FIG. 3: The qualitative form of the pseudo-potential V3 in Eq. 
g > 1 (on the right). 



for 5 < 1 (on the left) and 



of the purely classical pseudo-potential, the equation for the critical points read 

n = - s'-s^ + ^ = 0, (45) 

where 

480 7r3rc^ 

9 = j^2.A (46) 



4 

Fe 



is a new dimensionless parameter. In deriving Eq. (I45l) . it was omitted a term negligible 
except if s ~ c^/vp^, which is unlikely. 

The quantity g plays a decisive role on the shape of V3. Indeed, calculating the discrim- 
inant shows that the solutions to the cubic in Eq. fj45l) are as follows: {a) g < 1 —>■ three 
distinct real roots (one negative and two positive); (b) (7 = 1 — > one negative root, one 
(positive) double root; (c) g > 1 ^ one (negative) real root, two complex conjugate roots. 
Therefore, g < 1 is the condition for the existence of a potential well, which can support 
oscillations. This is shown in Figure 3. The analytic formulae for the solutions of the cubic 
in Eq. fl45l) are cumbersome and will be omitted. 

Restoring physical coordinates, the necessary condition for oscillations is rewritten as 



9<l - || |ENr>^met;^eC, (47) 

where 7 = e"^ /ineohc ~ 1/137 is the fine structure constant. From Eq. (1471) it is seen 
that for sufficient electrostatic energy the width cr of the localized envelope field can show 
oscillations, supported by the competition between classical refraction and quantum diffrac- 
tion. Also, due to the Fermi pressure, for large particle densities the inequality fH7|) becomes 
more difficult to be met, since vpe ~ '"'o''^- For example, when uq ~ 10^^ (white dwarf), 
the right-hand-side of Eq. (H7|) is 0.6 GeV. For uq ~ 10^^ (the next generation intense 
laser-solid density plasma experiments), it is 57.5 MeV. 
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Finally, notice that TCad,sc from Eq. (120|) . evaluated with the variational solution (l36l) . 
is proportional to + V3, which is a constant of motion for Eq. fH3|) . Therefore, the 
approximate solution preserves one of the basic first integrals of the vector NLS equation 
f ITS]) . as it should be. 

VI. CONCLUSION 

In this paper, the quantum Zakharov system in fully three-dimensional space has been 
derived. An associated Lagrangian structure was found, as well as the pertinent conservation 
laws. From the Lagrangian formalism, many possibilities are opened. Here, the variational 
description was used to analyze the behavior of localized envelope electric fields of Gaussian 
shape, in both two- and three-space dimensions. It was shown that the quantum correc- 
tions induce qualitative and quantitative changes, inhibiting singularities and allowing for 
oscillations of the width of the Langmuir envelope field. This new dynamics can be tested 
in experiments. In particular, the role of the parameter g and the inequality in Eq. (HTll 
should be investigated. However, the variational method was applied only for the adiabatic 
and semiclassical case, which allows to derive the quantum modified vector NLS equation 
( ITSjl . Other, more general, scenarios for the solutions of the fully three-dimensional quantum 
Zakharov system are also worth to study, with numerical and real experiments. 

Acknowledgments 

This work was partially supported by the Alexander von Humboldt Foundation. Fernando 
Haas also thanks Professors Mattias Marklund and Gert Brodin for their warm hospitality 
at the Department of Physics of Umea University, where part of this work was produced. 



[1] V. E. Zakharov, Zh. Eksp. Teor. Fiz. 62, 1745 (1972) [Sov. Phys. JETP 35, 908 (1972)]. 

[2] M. V. Goldman, Rev. Mod. Phys. 56, 709 (1984). 

[3] S. G. Thornhill and D. ter Haar, Phys. Reports 43, 43 (1978). 

[4] L. G. Garcia, F. Haas, L. P. L. de Oliveira and J. Goedert, Phys. Plasmas 12, 012302 (2005). 

[5] F. Haas, G. Manfredi and M. R. Feix, Phys. Rev. E 62, 2763 (2000). 



15 



[6] G. Manfredi and F. Haas, Phys. Rev. B 64, 075316 (2001). 

[7] F. Haas, Phys. Plasmas 12, 062117 (2005). 

[8] M. Marklund, Phys. Plasmas 12, 082110 (2005). 

[9] F. Haas, Phys. Plasmas 14, 042309 (2007). 

[10] P. K. Shukla and B. Eliasson, Phys. Rev. Lett. 96, 245001 (2006); Phys. Lett. A 372, 2893 
(2008). 

[11] X. Y. Tang and P. K. Shukla, Phys. Scripta 76, 665 (2007). 

[12] M. A. Abdou and E. M. Abulwafa, Z. Naturforsch. A 63, 646 (2008). 

[13] S. A. El-Wakil and M. A. Abdou, Nonl. Anal. TMA 68, 235 (2008). 



[14] Q. Yang, C. Q. Dai, X. Y. Wang and J, 
the comments about this work in Ref. 



'. Zhang, J. Phys. Soc. Japan 74, 2492 (2005). See 



[15] A. P. Misra, D. Ghosh and A. R. Chowdhury, Phys. Lett. A 372, 1469 (2008). 

[16] V. E. Zakharov, A. F. Mastryukov and V. H. Sinakh, Fiz. Plazmy 1, 614 (1975) [Sov. J. 

Plasma Phys. 1, 339 (1975)]. 
[17] V. E. Zakharov, Handbook of Plasma Physics, eds. M. N. Rosenbluth and R. Z. Sagdeev 

(Elsevier, New York, 1984), vol. 2, p. 81. 
[18] M. Landman, G. C. Papanicolaou, C. Sulem, P. L. Sulem and X. P. Wang, Phys. Rev. A 46, 

7869 (1992). 

[19] G. C. Papanicolaou, C. Sulem, P. L. Sulem and X. P. Wang, Phys. Fluids B 3, 969 (1991). 
[20] D. F. Dubois, A. Hanssen, H. A. Rose and D. Russel, J. Geophys. Res. 98, 17543 (1993). 
[21] P. A. Robinson and D. H. Newman, Phys. Fluids B 2, 3120 (1990). 
[22] P. Y. Cheung and A. Y. Wong, Phys. Fluids 18, 1538 (1985). 

[23] H. Alinejad, P. A. Robinson, I. H. Cairns, O. Skjaeraasen and C. Sobhanian, Phys. Plasmas 
14, 082304 (2007). 

[24] K. Akimoto, H. L. Rowland and K. Papadopoulos, Phys. Fluids 31, 2185 (1988). 

[25] L. H. Li and X. Q. Li, Phys. Fluids B 5, 3819 (1993). 

[26] G. Pelletier, H. Sol and E. Asseo, Phys. Rev. A 38, 2552 (1988). 

[27] L. Stenflo, Phys. Rev. Lett. 48, 1441 (1982). 

[28] M. Marklund, G. Brodin and L. Stenflo, Phys. Rev. Lett. 91, 163601 (2003). 

[29] J. Gibbons, S. G. ThornhiU, M. J. Wardrop and D. ter Haar, J. Plasma Phys. 17, 153 (1977). 

[30] B. Malomed, D. Anderson, M. Lisak, M. L. Quiroga-Teixeiro and L. Stenflo, Phys. Rev. E 55, 

16 



962 (1997). 

[31] R. Fedele, U. de Angelis and T. Katsouleas, Phys. Rev. A 33, 4412 (1986). 



17 



